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Abstract. 

Elements heavier than hydrogen or helium that are present in the atmospheres of white dwarfs 
with effective temperatures lower than 25,000 K, are believed to be the result of accretion. By mea- 
suring the abundances of these elements and by assuming a steady-state accretion, we can derive 
the composition of the accreted matter and infer its source. The presence of radiative levitation, 
however, may affect the determination of the accretion rate. We present time-dependent diffusion 
calculations that take into account radiative levitation and accretion. The calculations are performed 
on C, N, O, Ne, Na, Mg, Al, Si, S, Ar, and Ca in hydrogen-rich white dwarf models with effective 
temperatures lower than 25,000 K and a gravity of \ogg — 8.0. We show that in the presence of 
accretion, the abundance of an element supported by the radiative levitation is given by the equilib- 
rium between the radiative and gravitational accelerations, unless the abundance predicted by the 
steady-state accretion is much greater than the abundance supported by the radiative acceleration. 
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INTRODUCTION 

The accretion of circumstellar matter is one of the most plausible physical processes 
that can explain the appearance of metals observed in the atmospheres of several white 
dwarfs with effective temperatures lower than 25,000 K (see, e.g., Dupuis, Fontaine, & 
Wesemael [1] and Koester & Wilken [2]). In this scenario and under the assumption of 
a steady-state, the matter that falls onto the surface of a white dwarf diffuses downward 
because of the gravitational settling, and at the same time, it is continuously replenished 
with new material in such a way that the observed abundances remain constant. This 
accretion-diffusion model can maintain appreciable traces of metals in the atmospheres 
of white dwarfs and account for the observed abundances. The discovery of circumstellar 
dusty and gas disks, or both, around the most metal-rich white dwarfs goes in the 
direction that these disks are the source of the accreted matter (see, e.g., Zuckerman 
& Becklin [3]; Reach et al. [4]; Gansicke et al. [5]; Farihi, Jura, & Zuckerman [6]; and 
Brinkworth et al. [7]). Spectroscopic optical observations of such white dwarfs suggest 
that the circumstellar disks are the remaining of extrasolar minor planets that have been 
tidally disrupted by the white dwarf (see, e.g., Zuckerman et al. [8]; Klein et al. [9]; and 
Dufouretal [10]). 

The accretion-diffusion model shows that under the steady-state assumption, the 
abundances of metals observed in the atmospheres of cool white dwarfs can be related 



to the composition of circumstellar disks. For this reason, white dwarfs are excellent 
probe to determine the composition of debris disks, and ultimately, the composition 
of the extrasolar minor planets. The accretion-diffusion model is only valid, however, 
when no other physical mechanisms interfere with the accretion. In this paper, we 
demonstrate that the radiative levitation on elements such as carbon, aluminum and 
silicon in relatively cool DA models is important and cannot be neglected when one 
tries to determine the composition of the accreted matter. 



RADIATIVE LEVITATION 



Figure 1 shows the radiative acceleration as a function of depth for 1 1 elements in a DA 
model with T e ff = 20,000 K and logg = 8.0. These calculations are based on the work 
of Chayer et al. [11]. The solid curve gives the maximum radiative acceleration that 
is transferred to an element when all its absorption lines are desaturated. On the other 
hand, by increasing the abundance the saturation of the lines increases and, as a result, 
the radiative acceleration decreases. This is illustrated by the dotted curve in Figure 1. 
The radiative acceleration in this case is computed for an abundance of mass fraction 
X(A) = 10~ 2 . The figure shows that aluminum and silicon have radiative accelerations 
that are greater than the gravitational acceleration (dashed curve). Consequently, these 
two elements can be maintained in the atmosphere of such a white dwarf by the radiative 
levitation. 




FIGURE 1. Radiative acceleration as a function of depth for 1 1 elements in a DA model with T e s = 
20,000 K and logg = 8.0. The solid curve gives the maximum radiative acceleration that can be transmitted 
to an element, i.e., in the regime where the lines are completely desaturated. The dotted curve gives the 
radiative acceleration for an abundance of mass fraction X(A) = 10~ 2 . The horizontal dashed curve gives 
the local effective gravity. 
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FIGURE 2. Evolving distribution of Si in a DA model in the presence of a weak accretion, but without 
considering the effect of the radiative acceleration. Silicon is accreted at a rate of M acc = 100 g s _1 onto 
the surface of a DA model with T e ff = 20,000 K and logg = 8.0. Each profile is labeled by an increasing 
number that gives on the right hand side of the plot the corresponding time in years. We considered a 
uniform Si abundance of logiV(Si)/iV(H) = —7.0 as the initial condition. 



In order to answer the question on how the radiative levitation affects the calculations of 
accretion rates, we have carried out two sets of time-dependent diffusion calculations. 
The first set includes diffusion in presence of accretion alone, and the second set includes 
diffusion in presence of both accretion and radiative levitation. The diffusion of an 
element in an white dwarf envelope in presence of accretion is followed by resolving the 
continuity equation with a boundary condition that acts as a source term at the surface of 
the star. The radiative acceleration, on the other hand, is considered as an additional term 
in the microscopic diffusion velocity that enters in the continuity equation. Examples of 
time-dependent diffusion calculations can be found in Chayer, Fontaine, & Pelletier [12], 
Dupuis et al. [13]. and Vennes et al. [14]. 

Figure 2 shows the diffusion of Si as a function of depth in presence of accretion 
in a DA model with T e ff = 20,000 K and logg = 8.0. The calculations start with an 
homogeneous Si abundance corresponding to log./V(Si)/iV(H) = —7.0 (curve labeled 
by the number 1), and by considering a continuous accretion of Si onto the surface 
of the white dwarf at a rate of M acc = 100 g s . The top of the envelope is at the 
Rosseland optical depth Tr = 2/3. Although the calculations are performed for a period 
of t = 10 6 years, the figure shows the distribution of Si at 10 time steps from t = 
to t = 10,180 years. The figure shows that under the influence of the gravitational 



ACCRETION AND RADIATIVE LEVITATION 





2) 0.01 

3) 0.1 

4) 0.5 

5) 2.3 

6) 8.7 

7) 37. 

8) 160. 

9) 646. 



1) 0.00 yrs 



10) 2597. 

11) 10292. 



16 



14 



-12 

log AM/M 



10 



8 



6 



FIGURE 3. Same as Fig. 2, but the radiative acceleration on silicon is taken into account. 



settling, Si diffuses very rapidly from an abundance of logiV(Si)/iV(H) = —7.0 to an 
abundance corresponding to the steady-state value, i.e., the abundance for which there 
is an equilibrium between the accretion and the gravitational settling. The steady-state 
abundance at Tr = 2/3 is logiV(Si)/iV(H) = —10.8. At greater depths, as Si continues 
to flow downward it reaches a steady state abundance at each depth. 

Figure 3 shows the evolution of the Si distribution as we have discussed in the previous 
paragraph, but this time we consider the effect of the radiative levitation. Silicon is still 
accreted at rate of M acc = 100 g s _1 . We have seen in Figure 1 that Si is maintained 
in the upper layers of the atmosphere of a DA white dwarf with T e ff = 20,000 K 
and logg = 8.0. According to the radiative levitation theory, Si is supported because 
the upward radiative acceleration is equal to the downward gravitational acceleration 
(Chayer et al. [11]). By taking into account the radiative acceleration in the time- 
dependent diffusion calculations, we show in Figure 3 that Si diffuses rapidly and then 
reaches the equilibrium abundance given by the radiative levitation. A reservoir of Si 
builds up in the region where the radiative acceleration is greater than the gravitational 
acceleration. At deeper depths, however, where no radiative support is possible, Si flows 
downward and reaches a steady state that is almost identical to the one illustrated in 
Figure 2. The slight difference is due to the radiative acceleration that slows down the 
diffusion of Si. In contrast with Figure 2, the Si abundance on top of the envelope 
is no longer the abundance given by the steady state, but it is rather given by the 
radiative levitation. The Si abundance at Tr = 2/3 is iV(Si)/iV(H) = —8.1, which is 
about 500 times greater than the steady-state abundance when Si is accreted at a rate of 
M acc = 100 gs- 1 . 

We have computed the evolution of the Si distribution in a DA model when Si is 
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FIGURE 4. Comparison of the distributions of Si in a DA model in the presence of accretion computed 
at different accretion rates to models that include the radiative acceleration on Si (solid curve), and to 
models that do not include it (dotted curve). These distributions are obtained by following the diffusion of 
Si throughout the envelope of a DA model with r e ff = 20,000 K and \ogg = 8.0 for 10 6 years. The curves 
are labeled by the logarithmic value of the accretion rate M acc in g s _1 . 



accreted onto the surface of the star at a very low accretion rate, M acc = 100 g s _1 . What 
happens to the Si distribution when Si is accreted onto the star at a higher rate? How 
the equilibrium abundance supported by the radiative levitation is affected by higher 
accretion rates? Figure 4 illustrates this situation by comparing the distributions of Si 
in our DA model when Si is accreted at different accretion rates. The distributions of 
Si are computed by taking the radiative acceleration into consideration (solid curve), 
and without taking it into consideration (dotted curve). The accretion rates range from 
M acc = 10 2 to 10 7 g s -1 . Each curve in the figure is labeled by the logarithmic value of 
the accretion rate. The figure shows that the distributions of Si that are computed with 
and without the radiative levitation are almost identical for depths log AM/M —15. 
However, in the regions where the radiative support becomes important and when Si 
is accreted with rates of M acc ^ 10 4 g s _1 , the abundance at the top of the envelope is 
not set by the steady-state abundance, but by the equilibrium abundance supported by 
the radiative levitation. On the other hand, when Si is accreted at higher rates, i.e., for 
accretion rates M acc ~ 10 5 g s _1 , the Si abundances at the surface of the star converge 
to the steady-state abundances. For these high accretion rates, the abundances are no 
longer given solely by the radiative levitation, but rather by the equilibrium between the 
accretion and the downward diffusion. 



CONCLUSION 



Our findings demonstrate that the radiative levitation affects the observed abundances 
at the surface of DA white dwarfs when accretion is present. Although we have only 
showed time-dependent diffusion calculations on Si in one DA model (T e ff = 20,000 K 
and logg = 8.0), we have in fact carried out our investigation for a range of effective 
temperatures from 10,000 K to 25,000 K, one gravity logg = 8.0, and 1 1 elements. With 
this in mind, we can report that among the 1 1 elements that we have analyzed, C, Al 
and Si are the elements that are susceptible to the radiative support. If these elements 
are observed in the atmospheres of DA white dwarfs with 17,000 K £ 7kf ;$ 25,000 K, 
their abundances could be interpreted according to three scenarios: 1) when there is no 
accretion, the observed abundances should correspond to the radiative equilibrium abun- 
dances; 2) when there is accretion with low accretion rates, the observed abundances 
should still be the radiative equilibrium abundances; 3) when there is accretion with 
high accretion rates, the observed abundances should approach the steady-state abun- 
dances. The non-detection of these elements, however, could indicate that they have left 
the atmosphere of the white dwarf during the earlier phases of its evolution. 

For a direct application of the current model to a handful of DA white dwarfs with 
T e ff £ 25,000 K, see Dupuis, Chayer, & Henault-Brunet (these proceedings). 
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